Fermi-Bose mapping and N-particle ground state of spin-polarized fermions in tight 

atom waveguides 
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A K-matrix for wave-guide confined spin-polarized fermionic atoms recently computed by Granger 
and Blume is identified, in the low-energy domain, with a contact condition for one-dimensional (ID) 
spinless fermions. Difficulties in consistently formulating the contact conditions in terms of inter- 
action potentials are discussed and a rigorous alternative variational reformulation is constructed. 
A duality between ID fermions and bosons with zero-range interactions suggested by Cheon and 
Shigehara is shown to hold for the effective fD dynamics of a spin-polarized Fermi gas with 3D 
p-wave interactions and that of a Bose gas with 3D s-wave interactions in a tight waveguide. This 
generalizes the mapping from impenetrable bosons (TG gas) to free fermions and is used to derive 
the equation of state of an ultracold spin-polarized fermionic vapor in a tight waveguide. Near a ID 
confinement-induced resonance one has a "fermionic TG gas" which maps to an ideal Bose gas. 
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Ultracold atomic vapors in atom waveguides are cur- 
rently a subject of great experimental and theoretical in- 
terest and activity due to potential applicability to atom 
interferometry [1, 2] and integrated atom optics [3, 4] and 
their utility for demonstrating novel highly-correlated 
quantum states. Exploration of these systems is facili- 
tated by tunability of their interactions by external mag- 
netic fields via Feshbach resonances [5]. In fermionic 
atoms in the same spin state, s-wave scattering is forbid- 
den by the exclusion principle and p-wave interactions 
are usually negligible. However, they can be greatly en- 
hanced by Feshbach resonances, which have recently been 
observed in an ultracold atomic vapor of spin-polarized 
fermions [6]. Additional resonances are induced by tight 
transverse confinement in an atom waveguide. Of partic- 
ular interest is the regime of low temperatures and den- 
sities where transverse oscillator modes are frozen and 
the dynamics is described by an effective ID Hamilto- 
nian with zero-range interactions [7, 8] , a regime already 
reached experimentally [9-11]. Transverse modes are still 
virtually excited during collisions, leading to renormal- 
ization of the effective ID coupling constant gio via a 
confinement-induced resonance. This was first shown for 
bosons [7] and recently explained in terms of Feshbach 
resonances associated with bound states in closed, virtu- 
ally excited transverse oscillation channels [12]. Recently 
the analogous problem for fermions has been solved by 
Granger and Blume [13], who have shown that such reso- 
nances also occur in spin-polarized fermionic vapors. In- 
vestigation of such systems is facilitated by a mapping 
which allows reduction of strongly interacting fermions in 
one dimension to weakly interacting bosons. An energy- 
dependent mapping of this type was demonstrated in this 
recent work [13]. The analysis herein will be limited to a 
low-energy regime where one can use a simpler mapping 
originally employed to reduce the ID hard core Bose gas 



to an ideal Fermi gas [14, 15]. Here, following Cheon 
and Shigehara [16], the same mapping will be employed 
to map the strongly interacting Fermi gas to a weakly 
interacting Bose gas. More generally, for all values of 
the effective ID fermionic coupling constant, the known 
ground state of the ID Bose gas with delta-function re- 
pulsion, the Lieb-Liniger (LL) gas [17], will be mapped 
to generate the ID ground state of a spin-polarized Fermi 
gas with zero-range p-wave interactions. 

Contact condition for spin-polarized fermions in a 
waveguide: Granger and Blume derived the effective one- 
dimensional K-matrix for two interacting fermions con- 
fined in a single-mode harmonic atom waveguide [13]. It 
can be shown that in the low-energy [18] domain the K- 
matrix can be reproduced, with a relative error as small 
as 0{kl), by the contact condition 
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is the odd-wave one-dimensional scattering length, Vp = 
Op = — limfc^o tan(5p(/c)/A:^ is the p-wave "scatter- 
ing volume" [19], Op is the p-wave scattering length, 
a± = y/h/Ji(jj± is the transverse oscillator length [20], 
C(-l/2,l) = -C(3/2)/47r = -0.2079... is the Hurwitz 
zeta function evaluated at (—1/2,1) [21], and n is the 
reduced mass. The expression (2) has a resonance at a 
negative critical value Vp"*/aj_ = —0.4009 • • ■. In accor- 
dance with (1), the low-energy fermionic wavefunctions, 
Eq. (20) of [13], are discontinuous at contact, but left 
and right limits of their derivatives coincide. Following 
[16] we assume the same here. 

Odd-wave one- dimensional interaction potential: Fol- 
lowing the even-wave (bosonic) case, where the d- 



2 



interaction can be introduced naturally to cancel the S- 
functions resulting from double-differentiation of func- 
tions with discontinuous derivatives, in the case of 
fermions whose wave function is discontinuous it is 
tempting to introduce 6' interactions. However, 6'- 
functions and second derivatives are known to be ill- 
defined if used in a convolution with discontinuous ftmc- 
tions, making a consistent Hamiltonian formulation and 
corresponding perturbative treatments impossible. How- 
ever, a consistent variational formulation does exist, 
where matrix elements of operators are replaced by 
two-slot functional not factorizable as standard "bra- 
operator-ket" products. Such a formulation docs allow 
an accurate first order perturbation theory, and higher 
orders are under investigation. Furthermore, an exact 
Fermi-Bose mapping to be discussed allows nonpertur- 
bative treatment in the equivalent bosonic space. 
Consider general contact conditions 

^'(0+) - V>'(0-) = -(«f^)-i[^(0+) + V'(O-)] 

V(0+) - V(O-) = -af^[^'(0+) + V'(O-)] (3) 

that can scatter both even and odd partial waves. Here 
afjj (ofu) is the even (odd) scattering length. Our goal 
is to identify a functional whose extrema are solutions of 
the free-space Schrodinger equation subject to the con- 
tact conditions (3). Introduce "two-slot" Hermitian func- 
tionals corresponding to the "square-derivative" and S- 
function respectively: 

+i[x*(o+)-x*(o-)][^'(o+) + ^'(o-)] 



+2 [x*'(o+) + x*'(o-)][V(o+) - ^(0-)] 



(4) 



I dzx* 5{z) ^ = \[x*{0+) + x*(0-)][V(0+) + ^(0-)] 

After a lengthy but straightforward calculation one can 
show that extrema ijj of the energy functional 

£ = h'^/2ii j dzi)*'^' + gfo j dzii)*5{z)i) 

+g(j,Jdzr'S{z)i>' (5) 

with integrals defined by (4) with x = ^'^d the vari- 
ational space spanned by (normalized) wave functions 
with arbitrary discontinuities at zero, do obey the con- 
tact conditions (3), being local eigcnstatcs of the ki- 
netic energy outside of the contact point z = 0. Here 
the coupling constants are given by g^jj = —h^/jia^jj 
and g^jj = +fi^afjj/iJ,. One may introduce a formal 
"Hamiltonian" for the relative motion of two fermions by 
H[o = (ftV2/i)-(a,)(5,)^+.9f^ ^{d,) S{z) (5,)^ where 
this "operator" (and especially the kinetic energy part of 



it) must never appear outside of matrix elements, which 
should be carefully computed using the rules (4), and 
the eigenvalue problem for this Hamiltonian must be re- 
placed by a variational one. Notice that according to (4), 
the kinetic energy operator {h? /2^)^ {dz){dz)^ is a "reg- 
ularized kinetic energy" defined in such a way that the 
product of two 5 function contributions is automatically 
subtracted from the result of insertion of this operator be- 
tween two discontinuous functions. We have verified that 
for two fermions in an anti-periodic box the "potential" 
giD^(dz) 6{z) {dz)~* correctly reproduces the first order 
perturbation theory correction to the energy, but we warn 
the reader that the formal similarity between functionals 
(4) and matrix elements of real operators should not lead 
to an attempt to reformulate the problem as a matrix di- 
agonalization with some basis set. For example, one can 
check that in momentum space such a procedure leads to 
ultraviolet divergences, and an attempt to cancel them 
leads again to an unfactorizable functional. However, 
the contact conditions plus the free-particle Schrodinger 
equation for ^ 7^ do define a well-posed eigenvalue prob- 
lem not requiring use of a formal interaction operator. 

Fermi-Bose mapping: On the space of antisymmet- 
ric functions tl)p the contact conditions (3) reduce to 
^-^(0+) = -^f{Q-) = -afoVF(0±) with ^'p{Q+) = 
V'ir(0— ), and on the space of symmetric functions 
they reduce to -05(0+) = -V'b(O-) = -(«fD)~ Vb(0±) 
with ^Pb{^+) = V'b(0~)- Defining symmetric wave 
functions tjjB{z) = sgn{z)'tpF{z) and mapped scattering 



length a^jj 



''ID 



= aiD where sgn(z) is if z > 



and —1 if 2; < 0, one finds that the Bose and Fermi 
contact conditions are equivalent. Since the kinetic en- 
ergy contributions from z 7^ also agree, one has a 
mapping from the fermionic to bosonic problem which 
preserves energy eigenvalues and dynamics. The rela- 
tion between coupling constants g^j^ in H^jj and gfjy 
in iJf^ = -{n^/2,,)dl+g'^j,5{z) is gfj, = -h^ 1 1x^1^, 
and by (2) this agrees with the low-energy limit of Eq. 
(25) of [13, 18]. In the limit = +oc arising when 
0—, this is the N = 2 case of the original 
mapping [14, 15] from hard sphere bosons to an ideal 
Fermi gas, but now generalized to arbitrary coupling con- 
stants and used in the inverse direction. This general- 
izes to arbitrary N: Fermionic solutions ippizi, - " ^ ZN'-,t) 
are mapped to bosonic solutions iI)b{zi, ■ ■ ■ , zj^i^t) via 
tl^B = A{zi,- ■ ■ ,ZN)i>F{zi,- ■ ■ ,ZN;t)ipF where A = 
Yli<j<i<N ^S^izji) is the same mapping function used 
originally [14, 15]. The Fermi contact conditions 
are tpF\zj=ze-^- = -"^rlzj^ze- = -{aiD/2){dz^ - 
dze)4>F\zj=ze± and imply the Bose contact conditions 

{dz^ ' d^,)lpB\zi=z,+ = -{dz, - dz,)^JB\z,=z,- = 



-{2/aiD)i'B\zj=ze with aiD 



'■ID 



'-ID' 



and 



these are the usual LL contact conditions [17]. This 

mapping remains valid if external potentials Vext{zj) 
and/or additional interactions vi,r.{zji) of nonzero range 
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are present. One can define a formal fermionic 
Hamiltonian Hf^ = (?iV2m) Ef=i + 
9iD Y.i<j<£<N ^i9zji)S{z) {dzjt)^, but we again warn 
tiie reader that in calculations it should be treated varia- 
tionally or its "interaction" term replaced by the contact 
conditions, and it must not be substituted into a second- 
quantized framework. On the other hand, the definition 
of the interaction term in H^j-, = — (/i^/2/z) X^jLi + 

9?D J2i<j<£<N ^{zji) is much less delicate, so calcula- 
tions, including second quantization if desired, can be 
performed in the mapped bosonic Hilbert space. 

N-particle states: The exact ground [17] and excited 
[22] states of Hfjj are known for all positive gfjj if 
no external potential or nonzero range interactions are 
present, and the mapping then generates the exact N- 
body ground and excited states of H[jj. Define di- 
mensionless bosonic and fermionic coupling constants by 
jB = rngijj/nh^ and = — wgfpn/fi^ where n is the 
longitudinal particle number density and the minus pref- 
actor of 7i? is convenient since gfjj and have opposite 
signs. They satisfy 7b7f = 4. The ground state energy 
per particle e is related to a dimensionless function 6(7) 
available online [23] via e = (?i^/2m)n^e(7) where 7 is 
related to 7^ herein by 7 = 73 = 4:/^p. This is plotted 
as a function of jf in Fig. 1. If 7ir and jb are negative 
the Bose gas and mapped Fermi gas are unstable against 
collapse to an ultrahigh density droplet ( "bright soliton" ) 
in the absence of longitudinal trapping [24]. However, 
the derivation of an effective ID Hamiltonian from the 
3D one breaks down in the collapsed regime, as does 
neglect of three-particle and in fact multiparticle inter- 
atomic interactions. In the case of longitudinal trapping 
the gaseous regime is probably metastable if |7_b| is not 
too large (hence |7f| not too small), but this is beyond 
the scope of the present treatment. 

Fermionic TG gas: The mapping V^b = A^jjp was orig- 
inally introduced to map the strongly-interacting many- 
body problem of ID hard-sphere bosons of diameter d 
to the ideal Fermi gas [14, 15]. The simplest case d 
0+, the impenetrable point ID Bose gas, has recently 
elicited a great deal of theoretical [25 29] and experi- 
mental [10, 11] activity in the context of bosonic atomic 
vapors in tight atom wave guides, where it is now called 
the Tonks- Girardeau (TG) gas [30 36]. For bosons the 
TG regime is reached when gfj-, is large enough and/or 
the density n low enough that 73 1. A similar simpli- 
fication occurs in the fermionic case, where a fermionic 
TG regime is reached when 5^ is negative and large 
enough and/or n high enough that 7^ 1. The corre- 
sponding fermionic TG gas then maps to the ideal Bose 
gas since 757^ = 4. As an example, suppose that there 
is a longitudinal trap potential Vext{z) = ("^/2)'^forag'^^- 
Then in the fermionic TG limit the iV-boson ground 
state is Vb(2i, • • • , -zat) = TljLi "0(2:2) with uo{z) = 
7r"^/'^a,^yg^e~(^/"'°'"')^ with aio„g = y/h/jjuiong, and the 
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FIG. 1: Log-log plot of scaled ground state energy per par- 
ticle e = 2me/h^n^ versus dimensionless fermionic coupling 
constant 7f. 

corresponding fermionic TG ground state tpF = AipB 
has discontinuities at collisions zj = zi. Fig. 2 shows 
V'o^ and for AT = 3. The discontinuities in ■^q' are 
a consequence of idealization to a zero-range pseudopo- 
tential. For a potential of nonzero range tq -C Op they 
are rounded over a distance <C a^. As an illustration. 
Fig. 3 compares the two-particle ground state of the un- 
trapped fermionic TG gas with the solution when the 
zero-range interaction is replaced by a square well po- 
tential equal to — Vio when —zq < z < zq and zero when 
\z\ > Zq. (Note that the interaction term in H[jy is nega- 
tive definite in the regime of interest, where gfj^ < and 
9 ID > 0. ) The energy is taken as zero so the exterior solu- 
tion is sgn(z) = ±1; an interior solution fitting smoothly 
onto this is sin(Kz) with k = a/2/xVo/^^ = 7r/2zo, the 
critical value where the last bound state passes into the 
continuum, a zero-energy resonance. A fermionic con- 
tact condition with a finite scattering length can be ob- 
tained in the limit 2:0 ^ if k scales with the width zq 
asK= {tt/2zo)[1 + (2/7r)2(zo/af^)]. 

Discussion: The effective ID A'-particle ground state 
of a spin-polarized Fermi gas with zero-range p-wave in- 
teractions has been mapped to the A^-particle ground 
state of the ID Bose gas with delta-function repulsion 
[17], providing the exact solution of this fermionic prob- 
lem in the absence of longitudinal trapping. Experiments 
on spin-polarized Fermi gases in this quasi-lD regime are 
suggested, as is investigation of Fermi-Bose duality for 
waveguide-confined ultracold gases with realistic interac- 
tions. 
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PHY-0301052 (M.O.). 
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FIG. 2: tpoFiz, Z2, Z[i) (solid lino) and ^ob{z, Z2, z;i) (dashed 
line) for a longitudinally trapped fermionic TG g func- 
tion of z for Z2 = 0.5 and za = 1- Units are such that 
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FIG. 3: iV = 2 untrapped fermionic TG gas ground state 
(dashed line) compared with zero-energy scattering solution 
for a square well with range zo and depth Vb corresponding to 
the boundary between no bound state and one bound state, 
a zero energy resonance (solid line), as function of relative 
coordinate z. Units are such that zo = 1. 
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